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Abstract. We present algorithms to factorize homogeneous polynomi- 
als in the first g-Weyl algebra and the first Weyl algebra. By homoge- 
neous we mean homogeneous with respect to a Z-grading on the first 
((j'-)Weyl algebra which will be introduced in this article. It turns out 
that the factorization can be almost completely reduced to commutative 
univariate factorization with some additional easy combinatorial steps. 
We im plemented the algor it hms in the computer alg ebra system Sin- 
gular (jPecker et all (|2012l '). [Greuel and PfisteJ (|2007l )'). which provides 
algorithms and an i nterface to deal with noncommutative polynomials 
l|Greuel et al.l l|2010l) ). The implementation beats currently available im- 
plementations dealing with factorization in the first Weyl algebra in 
speed and elegancy of the results. Furthermore it broadens the range 
of polynomials we can nowadays factorize in a feasible amount of time 
using a computer algebra system. 



1 Introduction 

Operator algebras as the g-Weyl and the Weyl algebras are interesting objects 
to study as one can derive properties of the solution spaces of their associated 
equations one wants to solve. Especially concerning the problem of finding the 
solution, the preconditioning step of factorizing an operator can help to reduce 
enormously the complexity of that problem in advance. 

But usually those operator algebras are noncommutative polynomial alge- 
bras, and a factorization of an element in those algebras is neither unique in the 
classical sense (i.e. unique up to multiplication by a unit), nor easy to compute 
at all in general. 

Nevertheless, a lot has been done in this field in the past. S. Tsarev has 
studie d the for r n, num ber and the prope rties of the factors of a differential oper- 
ator in Tsarev (Il994 ) and Tsarev (1 1996 ), wh ere he used and extended the work 
presented in iLoewvl ( 19031 ) and iLoewvl ( 19061 ). 

A very general approach to noncommutative algebras and their pr o pertie s, 
including factorization, is also done by J. Bueso et al. in (jBueso et al.l (|2003l) '). 
They provide several algorithms there and introduce various point of views when 
dealing with noncommutative polynomial algebras. 

In his dissertat ion M. van Hoeij developed an algorithm to factorize a dif- 
ferential operator (|van Hoeiil ^9M))- There were several papers fo llowing that 
dissertation using and extending those techniques (e.g. Ivan Hoeii 1 (Hmz,) and 



van Hoeii and Yuaiil ( 201Cll )). and nowadays this algo rithm is implemented in the 



DETools package of Maple (jMonagan et alj (|2Q08f )) as the standard algorithm 



for factorization of those operators. 

From a more algebraic point of view and dealing only with strict polyno- 
mial noncommutative algebras, H. Melen k and J. Apel developed a package for 
the computer algebra system REDUCE (iMelen k and Apell (ll994|)). That pack- 



age provides tools to deal with noncommutative polynomial algebras and also 
contains a factorization algorithm for the supported algebras. 

Those algorithms and implementations are very well written and they were 
able to factorize a large amount of polynomials we gave them as input. Nonethe- 
less, as we will see in this paper, there exists a large class of polynomials that 
seem to form the worst case for those algorithms. There, one can use a different 
approach t o obtai n a fac torization very quickly. This approach was originally de- 
veloped in iHeinle and further extended and presented at the conference 



ISSAC in Munich, 2010, the DEAM2 workshop in Linz, 2011, and the conference 
"Computeralgebra Tagung" in Kassel, 2012. Now, this is the first journal version 
of this work. Furthermore, the algorithms were implemented in Singular, and 
since version 3-1-3 they became part of the distribution (library; ncfactor.lib). 

1.1 Basic Notions and Definitions 

We will start with introducing the first q-Weyl algebra and the first Weyl algebra. 
By K we always denote an arbitrary field. 

Definition 1.1. The first q-Weyl algebra Qi is defined as 

Qi := K{x,d\dx = qxd+l), 

where q is a unit in K. Qi is the operator algebra associated to 

f{q^) - fix) 



dq ■■ fix) ^ 



{q-l)x 



also known as the q-de rivative, where f G K.[x] . For further reading consider 
Kac and Cheund l200A ). 



For q = \, the operator is still well defined. This can he seen in the following 
way. Let f = X]r=o '^iX^ j where n G Nq and a; G K. Then 

n n n 

fiqx) - fix) ^ ^a,{qxy - ^a^ = ^aix^q"- - 1). 

i=0 i=0 1=0 

The expression q — 1 is clearly a divisor of q^ ~ I for all i >1, and we obtain 



fjqx) - fix) ^ i-i I 

iq-l)x ^Z."'" [l-^ 

^ i=l \j=Q 



For the special case where q = I we have the first Weyl algebra, which is 
denoted by Ai . 
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The first (g-) Weyl algebra possesses a. nontrivial Z-grad i ng - introduced b y 
M. Kashiwara and B. Malgrange in 1983 ()Kashiwaral (|l983t ). iMalgrangd (|l983l )') 
- using the weight vector [— for a v G Z on the tuple [x, i9]. For simplicity, 
we will choose v := 1. In what follows, deg denotes the degree induced by this 
weight vector. By homogeneous polynomials, we also mean homogeneous with 
respect to that weight vector. 

Example 1.2. We have 

deg(a;9) = deg(9a; + 1) = 0. 

The polynomial 

xd'^ + x^d"^ + d 
is also a homogeneous polynomial, but of degree one. 

The nth graded part of Qi and analogously the nth graded part of Ai is 
given by 



i.e. the degree of a monomial is determined by the difference of its powers in 
x and d. 

Concerning this choice of degree, the so called Euler operator 

e -.^xd 

will play an important role as we will see soon. 

First of all, let us study some commutation rules the Euler operator has with 
X and d. For Qi, in order abbreviate the size of our formulas, we introduce the 
so called g-bracket. 

Definition 1.3. For n £ N, we define the q-bracket [n]q by 



1 -g" 



J 9 



n-1 

1 - (7 E 



q 



Lemma 1.4 (Compare with Saito et al. ( 2000h ^. InAi, the following 
mutation rules do hold for n G N; 

fe" =:=x"(0 + n) 

6»9" = a"(e'-n). 

More general, in Qi the following commutation rules do hold for n G N; 

eon = ^ _ g'"+' - 1 
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Those rules follow via induction on n e N and we recommend to perform this 
induction in order to get experience with calculating in the first q'-Weyl algebra. 



Remark 1.5. With the help of the Lemma above one can also easily see that the 
so called shift algebra 

K(n, s\sn = (n + l)s) 

is a subalgebra of the first Weyl algebra Ai. The embedding of a polynomial 

p = X]r=oP«('^)'*" algebra, where pi e K[n], into the first Weyl algebra 

is done via the following map: 

n n 

i : K{n,s\sn = (n + l)s) ^ Ai, ^pi(n)s" ^ ^p,{e)d'^ . 

Therefore, the factorization techniques developed here can also be applied to 
the first shift algebra. 



Those commutation rules can of course be extended to arbitrary polynomials 
in 61. 



Corollary 1.6. Consider f{0) := / G K[0],e := xd. Then, in Qi, for allnGN 
we have 

/(0)x"=x"/(g"(^+[n],) 



whereas in Ai we have 



q \ g""^ 1—9 



f{e)d^ = a"/ {e-n), 

Those are the basic tools we need to explain our approach for factoring 
homogeneous polynomials in the first Weyl and the first g-Weyl algebra. 



2 A new Approach for Factoring Homogeneous 
Polynomials in the First (q-) Weyl Algebra 

The main idea of our factorization technique lies in the reduction to a commuta- 
tive univariate polynomial subring of Ai rcsp. Qi, namely K[0]. It appears that 
this subring is quite huge in the sense of reducibility of its elements in Ai resp. 

Qi. 
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2.1 Factoring homogeneous polynomials of degree zero 

The following lemma shows, that we can rewrite every homogeneous polynomial 
of degree zero in Ai and Qi as polynomial in K[9]. 

Lemma 2.1 (Compare with ISaito et aiT ( 2000h . Lemma 1.3.1). In Ai, 

we have the following identity for n eN: 

n-l 

x"5" = [|(6'-i)- 
In Qi, one can rewrite as element in K[0] and it is equal to 




n — 1 



i=0 

where denotes the ith triangular number, i.e. 



for all i e No . 

Therefore the factorization of homogeneous polynomials of degree zero can 
be done by rewriting the polynomial as element in K[9] and apply a commutative 
factorization on the polynomial, which is well implemented in every computer 
algebra system. 

Of course, this would not be a complete factorization, as there are still ele- 
ments irreducible in K[9], but reducible in Qi resp. Ai. An obvious example is 
9 itself. 

But fortunately there are only two monic polynomials irreducible in K[9], 
but reducible in Ai resp. Qi. This is shown by the following Lemma. 

Lemma 2.2. The polynomials 9 and ^ + ^ a'^e the only irreducible monic ele- 
ments in K.[9] that are reducible in Qi- For Ai, the polynomials 9 and 9 + 1 are 
the only irreducible monic elements in K\9] that are reducible in Ai. 



Proof. We will only consider the proof for Qi, as the proof for Ai is done in an 
analogue way. Let / G 'K[9] be a monic polynomial. Assume that it is irreducible 
in IK[0], but reducible in Qi. Let (/?, tp be elements in Qi with ipip = f. Then ip 
and -tp are homogeneous and e Q[~^\i^ e Q[^^ for a fc e Z \ {0}. Without 
loss of generality let k be positive. For k being negative we can use a similar 
argument. 
Then 

ip = ip{9)x'' 
^ = ij{9)d'' 
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for i^jip 

Using Corollary II .61 we have 



, 9 V 9" ^ 1 — g 
As we know from Lemma 11.41 the equation 

i 



^ A;-l i-l 



i=o ' j=a 



holds. 

Thus, because we assumed / to be irreducible in M.[9], we must have ip,ip gK. 
and k = 1. Because / is monic, we must also have ip — ijj^^ . 

As a result, the only possible / is / = 0. If we originally would have chosen 
k to be negative, the only possibility for / would be / = + i. This completes 
the proof. 

Therefore, we have an algorithm to factor a homogeneous polynomial p £ Ai 
(resp. p £ Qi) of degree zero in K[d]. It is done using the following steps. 

1. Rewrite p as an element in M.[9]. 

2. Factorize this resulting element in K[d] with commutative methods. 

3. If there is or 6* + 1 among the factors, factorize it to a; • 9 resp. d ■ x. 

As mentioned before, the factorization of a polynomial in a noncommutative 
ring is in general not unique up to multiplication by units or interchanges. Sev- 
eral factorizations can occur. Fortunately, in the case of the polynomial first (q-) 
Weyl algeb ra, ther e are o nly finitely many different factorizations possible due 



to Tsarev ( TsarevI (19961)). In order to obtain all these different factorizations. 



one can apply the commutation rules for x and d with 9. That these are all pos- 
sible factorizations up to multiplication by units can be seen using an analogue 
approach as in the proof of Lemma 12.21 Let us consider an example for that. 



Example 2.3. Let 

p ■= x^d^ + 4x^5^ -I- ixd e Ai 

and K := Q. It is homogeneous of degree zero and rewritten in IK[6'] it is equal 
to 

6*3 + 6*2 + 9. 

This polynomial factorizes in Wifi] to 9 ■ {9^ +6* + 1), which further factorizes as 9 
is a factor to a;- 9 - (6*^ + 6*4- 1). To get all possible factorizations of p, we apply the 
commutation rules with x resp. d and obtain the following other factorizations. 

(6*2 + 6* + 1) • X • a, 

X- {9^ +39 + 3)- d. 
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Remark 8.4- Let us consider the complexity of the first two steps mentioned 
above. 

Ad step 1: The polynomial p has due to the assumption of being homoge- 
neous of degree zero the form 



In order two transform it into an element in K[0], we have to apply the rewriting 
rule stated in Lemma |2. II for every term x^d^ in p. For that, one can make use 
of the fact that 

^n+lQU+l ^^ngn . (^Q _ ^ 

Thus, in order to perform step 1, we need to perform for every i G n a 
multiplication of a polynomial in K.[9] of degree i with a polynomial of degree 1. 

Ad step 2: Unfortunately, the factorization problem even in the univari- 
ate case does not have polynomial complexity. In general we face exponential 
complexity. But there are algorithms developed since the 1980s that appeared 
to be very fast in practice. One example of it is th e famous LLL a lgorit hm by 
Lenstra, Lenstra, Lavsz, Lszlo developed in 1982 (iLenstra et al. For 



further reading s on the complexity of the factorization problem we recommend 



KaltofenI ()1982[ ) 



2.2 Factoring homogeneous polynomials of arbitrary degree 

Fortunately, the hard work is already done and the factoring of homogeneous 
polynomials of arbitrary degree is just a small further step. 
The following well known Theorem reveals why. 

Theorem 2.5. Qf^ resp. A^f^ is a ring, generated by the element 9 :— xd as 
a ^-algebra. The other graded direct summands Q'^^^ resp. A''^^ are cyclic Q^^^ 
resp. A^^^ modules generated by the element x^^ , if k <0, or by , ifk>0. 

Proof. The first statement can be seen using Lemma [2Tl as we can identify Q^^^ 
resp. Af '^ with K[9]. 

(k) A (k) 

For the second statement recall that for a polynomial p E Q\ resp. p E A\' 
being homogeneous of degree k E li means, that the distance of the powers of x 
and d in every monomial of p equals k resp. —k. That means, we can divide by 
x~^ resp. from the right and obtain the desired module structure. 

Therefore, one factorization of a homogeneous polynomial p E Q\ resp. 

(k\ 

p E A\ oi degree fc G Z can be obtained using the following steps. 

1. Divide p by x~^ resp. 9*^ from the right. 

2. Factorize the remaining polynomial - which is homogeneous of degree zero 
- using the steps shown in the previous subsection. 
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This will lead to one factorization. Again, to obtain all possible factorizations, 
one simply swaps the factors homogeneous of degree zero with the ds resp. a;s. 
How this is precisely done is described in Algorithm 2. 

Now we can finally state the algorithm to factorize homogeneous polynomials 
in the first {q-)Wey\ algebra. 

To obtain one factorization, we simply use the Algorithm 1. 



Algorithm 1 HomogFac: Factorization of a homogeneous polynomial in the 
first (q-) Weyl algebra 

Input: h G A^™' (resp. h G Qi™"*), where m G Z 

Output: (/i, ...,/„) G yl? resp. (/i, . . . , /n) e Q?, such that /i ■...•/„ = ft, n G N 

1: if m 7^ then 
2: if m < then 

3: Get h G A^"' such that h = hx'"" 

4: factor := {x, . . . , x) 

— m times 

5: else 

6: Get h such that h = hd"" 

7: factor := (d, . . . , d) 

m times 

8: end if 

9: else 

10: h:=h 

11: factor := 1 

12: end if 

13: {f\, . . . , fi) := Factorization of h as element in ]K[0] {I G N) 

14: (/i, . . . , h := Substitute 61 by a; • 9 in (A, ... , /„) 

15: result :— () 

16: for i from 1 to Z do 

17: if fi=x-d then 

18: Append x and d to result 

19: else 

20: if fi = d x then 

21: Append d and x to result 

22: else 

23: Append fi to result 

24: end if 
25: end if 
26: end for 

27: Append each element in factor to result 
28: return result 



To obtain all factorizations, wc extend Algorithm 1 by Algorithm 2. 
The termination of those algorithms is clear, as we always only iterate over 
finite sets. The correctness follows by our preliminary work. 
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Algorithm 2 HomogFacAU: All factorizations of a homogeneous polynomial in 
the first (g-) Weyl algebra 

Input: h G A^™' (resp. h £ Qi™'), where m £ Z 
Output: {(/i, . ..,/„) e Ari/i •...•/„ = /i, n e N} 

1: {h,---,fy,g,---,g) ■■= HomogFac(^) 

{v £ No,g £ {a;, 9}, /i £ ^j"' concatenated, if needed} 
2: Rewrite each fi as element in K[0] 

3: result := {Permutations of (/i, . . . , f,^, g, . . . , g) with respect to the commutation 
rules} 

4: for . . . , (7„) £ result do 

5: for i from 1 to n do 

6: if gi = 6 (resp. gi — 9 + 1) then 

7: gi := X, d 

8: leftpart ■- {gi,.. . ,c/i_i,a;) 

9: rightpart := (9, gi+i, . . . , g„) 

10: swap s down in the first list and d up in the second 

11: put all possibilities of the left and the right part together and add them to 

result 



12 
13 
14 
15 



end if 
end for 
end for 
return result 



Remark 2.6. We made the general assumption that we are dealing with the 
polynomial first Weyl algebra. That means that we assumed the variables x and 
d as polynomial indeterminates. In practice one often finds functions that are 
given polynomial in 9, but rational in x. The factorization algorithm in Maple 
factorizes exactly those functions. They have some odd properties, as ther e might 
be infi nitely many distinct factorizations. But fortunately, as we studied in lHeinlj 
one can lift factorizations of a polynomial from the rational Weyl algebra 
to a polynomial in the Weyl algebra as defined in this paper. Furthermore, we 
can even say that if an element is irreducible in the polynomial firs t Weyl algebra , 
then it is also irreducible in the rational first Weyl algebra (see Heinle! ( 2012f ). 
Chapter 1, Lemma 3.13). 



3 Implementation and benchmarking 



We implemented the presented algorithms in Singular, and since version 3-1-3 
they are part of the distribution of Singular. The following example shows how 
to use the library containing them. 
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Example 3.1. Let h € Qi be the polynomial 

h := g^S^^lO^lO + ^16(^4 + ^3 ^ ^2 ^ ^ ^ ^^2^9q9 

+ 3g^2 + + Uq^° + 20q^ + 26q^ 
+30g^ + Slq^ + 26q^ + 20g^ + 13^^ + 7q'^ + 3q+ l)x^d^ 
^q\q^ + 2g« + + Qq^ + 7q^ + Sq^ + Qq^ + V + 2? + 1) 
(g'' + g-V + g + + g + 1)2:^9^ 
+g(g2 + g + + 2q'^ + 2g3 + 3g2 + 2g + 1) 
(g4 + g3 + g2 + 5 + l)(g2 + i)(^ + i)x^d^ 
+(g^° + 5g'' + 12g''^ + 21g' + 29g" + 33g^ 
+31g* + 24g^ + 15g^ + 7g + I2)x^d^ + Qx^d^ + 24 

and K = Q. We can use Singular to obtain all of its factorizations in the 
following way. 

LIB "ncf actor. lib"; 
ring R = (0,q) , (x,d) ,dp; 
def r = nc_algebra (q,l); 

setring(r) ; 

poly h = ... //See the polynomial defined above. 
hoinogfacFirstQWeyl_all(h) ; 
[1] : 
[1] : 
1 

[2] : 

x5d5+x3d3+4 
[3] : 

x5d5+6 

[2] : 
[1] : 
1 

[2] : 

x5d5+6 
[3] : 

x5d5+x3d3+4 

As one can see here, the output is a list containing lists containing elements 
in Q\. Those elements in Q\ are factors of h, and each list represents one possible 
factorization of h. 

If the user is interested in just one factorization the command homogf acFirstQWe 
instead of homogf acFirstQWeyl^ll can be used. The output will then be just 
one list containing elements in Qi. 

On the computer we used - 2 GB RAM, 2.33GHz Dual Core processor - this 
calculation needs 2.8 seconds. Compared to the factorization of 

{x^d^ + &){x^d^ + x^d^ + A) 
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as element in Ai, which takes less than a second, this seems to be way more slow 
considering that both algorithms have the same complexity. But this slowdown 
is not due to more steps that need to be done in the algorithm for the g-Weyl 
algebra, but due to the parameter q and the speed of calculating in Q{q) as the 
basefield instead of just in Q. 

In fact, there is no computer algebra system known to the authors that can 
factor polynomials in the first g-Weyl algebra Qi. Therefore, we cannot compare 
our algorithms in this case to other implementations. 

For the first Weyl algebra Ai, there are other implementations. We can draw 
a comparison to the DFactor method in the DETools package of Maple and the 
nc_f actorize_all method in the NCPoly library of REDUCE. We used version 
16 of Maple and version 3.8 of REDUCE. 

Remark 3.2. As mentioned before, the algorithm implemented in Maple fac- 
torizes over the rational Weyl algebra, i.e. the variable a; is a rational argument 
having adjusted commutation rules with d. This is a weaker assumption on the 
input since the ring that is dealt with there is larger. The comparison is still 
valid, as we are interested in the fact if a factorization can be found at all and - 
if so - whether we can lift it to a polynomial factorization without multiplying 
by xs. 

We will not go into detail about how the algorithm in Maple works, as it 
wo uld take s e veral pages to describe it. The interested reader can find details 
in Ivan Hoeiil (|l997l ). It works with collections of exponential parts and their 
multiplicities at all singularities of a given differential operator / and subsequent 
calculation of left and right hand factors. 

The algorithm implemented in REDUCE is also working with the polyno- 
mial Weyl algebra. In fact, the algorithm written there can be applied to a large 
amount of polynomial noncommutative rings. 

Details about the functionality of the algorithm in REDUCE are unfortu- 
nately not available. One can only try to understand it from the code that is 
given open source. It uses several Grobner Basis computations in order to find 
its solutions. 

Example 3.3. Consider again the element 

h := (x^a^ + 6) • {x^d^ + x^d^ + 4). 

— Singular: Found two factorizations in less than a second. 

— Maple: Found one factorization after 29 seconds; The factors are huge (size 
of the output file is around 100KB). 

— REDUCE: Did not terminate after 9 hours of calculation. 

Example 3.4. We experimented with other examples of homogeneous polynomi- 
als in the first Weyl algebra. The results are listed in the next table. An entry 
labeled with "- NT -" stands for "no termination after four hours", and the 
abbreviation "fcts" stands for "factorizations" . 
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Singular 



Maple 



REDUCE 



(x^'^a^'J + 5xd + 7)-x^- {x''d'' + 3x'd' +xd + 4): 



0.37s; 12 fcts. 



NT ■ 



NT 



TTT" 



{.r"0" + ()) ■ {,v"()-' + rUr + ■ 



4.8s; 132 fcts. |31.62s; 1 fcts. |- NT - 

(5a;'"a^" + Tx^'d'' + Sx'^d'' + 9x'd' + ex^g" + 5x^8^ + 8x^9^ + 5x^9^ + Qa;^^^ + 9xd + 6) • g^^: 



■NT 



0. 99s; 21 fcts. |- NT , 

(7x^^a^^ + x^^S^^ - x^'^d''^ - 3a;^"a^" + 2x«ay + + aj-'a'^ - a;^a^ - Qa;*^* + ^9 - 1)- 
(8x13513 + 3x1^^1^ + x"5" - 2xi°ai° + lOx^a^ - 3x^5^ + 2x^5^ + x-^a^ + 38x8 + 1) ■ d^: 

7tZ I ^T~ ' l^TTi 



36.02s; 504 fcts. 



NT 



(jjiugiu + 23x^a« + 3x«a« - 9x^9' - x^a*^ + 3x^5* + 6x^9^ + 4x5 + 1)- 
{-x^8^ + 4x^9^ - x^8^ + 4x^5^ - 5x^5^ + x^d"^ - 7x8 - 10) • x 



10. 



5.63s; 132 fcts. 



NT 



NT 



7^ 



(=25^ 



5^16516 _ 7.^15515 ^ 4^14^14 _ ^13gl3 ^ 3.12^12 _ 2x11511 



x^5^ 



5x^58 



6x^5^ + x-i^^ + 2x^8^ + 219.T252 + xd - 1) • (- 



552 



1524 - 32x235^3 ^ ^22^22^ 



7x^15^1 + 61.T^"5^" - 2.t1«51« 



,,.lf)fll6 



5"' + 2.t1-V91'^ - 2.T14514- 



125.43s; 230 fcts. 



NT 



NT 



(xi^Si^ + 13x^5^ - x«5« + 4x^5' + 13x^5" - 3x'^5'^ - 37x''5'' - x^S^ + x^5^ + x5 - 1)- 
(_^iogio _ 23x9 59 + 2,x^8^ + x^a^ - x^a^ - 2x^8^ - 2x^8^ + 2x^8^ - x^S^ - 2x5 - 2): 



0.27s: (i fctri 



XT 



NT 



(98xi^5i-'' + 40x115" + 98.X-13513 + 44xi^5i^ + 55xii5ii + 96xi"5i" + 95./;'^5'^+ 
7x^5^ + 56x^5^ + 50.Tf'5^' + 40x''^55 + 11x^5^^ + 40x^53 + 78x^5^+ 
13x5+19) • (61xi'^5i-'^ + 50x1^5" + 83x135" + llxi^Si^ + 89x^511+ 
55xi°5i" + 81.'j;°5'^ + 63x^5^ + 22x^5^ + 10x^^8^+ 
35x^5^ + 90x^5^ + 60x353 + 20x^5^ + 30x5 + 43): 



0.23s; 2 fcts. 



-NT ■ 



-NT 



(85x""5^'^ + 80xi''5i'' + 27xi«5i« + 74xi'5i' + 49x^'^5i'^ + 95xi^5i=' + 96xi''5i'* 
+37x13513 + 26x12512 + 93x11511 + 39xi°5io + 19x^59 + 48x^58 + 82x^5^ 

+26x^5*^ + 26x^5^ + 7x^5^ + 61x353 + 8x^5^ + 81x5 + 88)^: 

0.22s; 1 fcts. I- NT- I- NT - 



The conclusion we can draw at this point is: Even if homogeneous polynomials 

seem to be easy objects to factorizc, they seem to form a worst case class for the 
implementations in REDUCE and Maple. 

Therefore, with our algorithm we are now able to factorize more polynomials 
using computer algebra systems: In general homogeneous polynomials in Qi, 
and for Ai we have broaden the range of polynomials that can be factorized in 
a feasible amount of time or even sometimes at all. 
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4 Conclusion and Future Work 



With this paper, we contributed an algorithm and an implementation for fac- 
torization of homogeneous polynomials in the first q-Weyl algebra. 

Furthermore, we also considered the special case of the first Weyl algebra 
and showed that our algorithm beats for the large class of [—1, 1] homogeneous 
polynomials current implementations in terms of speed and elegance of the so- 
lutions. 

We can also construct a family of polynomials where the implementation in 
Singular is the only one that is able to factorize those elements in a feasible 
amount of time and memory consumption. 

The next thing to deal with would be gene ral polynomia ls in the first (g-) 
Weyl algebra. A first attempt to that we did in iHeinl We made highly 

use of our knowledge about the grading of the first Weyl algebra. The approach 
has been almost completely of combinatorial nature. Its speed and quality of 
solutions was comparable in many cases to the other implementations, but it 
was not elegant. It is also distributed with Singular since version 3-1-3. 

Recently, we developed inter alia some new tec hniques to factorize general 
polynomials in the first Weyl algebra in Heinle ( 20121) . They are way more elegant 
and provide a good performance and complexity. The implementation of those 
techniques will soon be committed to Singular and are distributed since version 
3-1-6. 

But still, there is a lot of work to do in that field and as we have seen in this 
paper, it is worth to consider special cases within the factorization problem as 
there are often way more elegant ways to deal with them. In our case, it was a 
very large class of polynomials. 

Another interesting question would be factorization in the nth Weyl algebra 
and how similar approaches can be applied there. 
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